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In a previous paper [Phys. Rev. D 83, 113013 (2011)] we have shown that the solar sector data 
(solar and KamLAND) are sensitive to the parameter #14, encoding the admixture of the electron 
neutrino with a fourth (essentially) sterile mass eigenstate. In that work we evidenced that such 
data prefer a non-zero value of #14 and that such a preference is completely degenerate with that 
of non-zero #13. In this report we show how the evidence of 813 > 0, recently emerged from global 
neutrino data analyses, lifts such a degeneracy and disfavors the case of sterile neutrino mixing. By 
excluding from our analysis the total rate information coming from the reactor experiments we untie 
our results from any assumption on their flux normalization. In this way, we establish the robust 
upper bound sin 2 8 14 < 0.04 at the 90% C.L. 

PACS numbers: 14.60.Pq, 14.60.St 



In a recent paper [l| we have introduced the theoretical 
framework needed to describe solar neutrino oscillations 
within the so-called 3 + s schemes endowed with s new 
sterile neutrinos (see also Q). In the same work, we 
have considered the constraints attainable within such 
schemes from the "solar sector" (solar and KamLAND 
data) showing that this dataset, while preferring a non- 
null admixture of the electron-neutrino with mass eigen- 
states "far" from the solar (^1,^2) doublet, is currently 
unable to distinguish if such a mixing is realized with 
the third standard mass eigenstate V3 or with new ones 
(1^3+1, ^3+s)- In the simplest 3 + 1 framework, this 
ambiguity translates into a degeneracy of the estimates 
of the standard mixing angle #13 and the new angle #14 
(see [l[ for the details of the parameterization of the lep- 
ton mixing matrix). 

After publication of [l| new data have been released 
that are relevant to the analysis therein performed. In 
particular, the long-baseline (LBL) accelerator experi- 
ments Tokai-to-Kamioka (T2K) 0] and the Main Injector 
Neutrino Oscillation Search (MINOS) [H have both evi- 
denced a phenomenon of — > v e conversion. Moreover, 
the reactor experiment Double-CHOOZ (D-CHOOZ) [|, 
currently operating only with the far detector, has found 
an indication of v e v e disappearance. These findings, 
if interpreted within the standard 3-flavor framework, 
point towards a non-zero value of #13, in line with the 
first indications arising from global neutrino data anal- 
ysis @ (see also @, [I]). In fact, with the inclusion of 
the new crucial piece of information, an updated global 
neutrino data analysis [i[ (see also [HI Ell) provides 1 
evidence of #13 > at more than 3c. 



1 The analyses in [si lld(l do not incorporate the D-CHOOZ result, 
whose inclusion would further reinforce the evidence of non-zero 
013 therein established. A preference for 613 > around the 3cr 
level has been shown also in the analysis performed in , which 
includes D-CHOOZ (together with MINOS and T2K), but not 
the solar and atmospheric data. 



This new circumstance prompts us to improve the 
analysis performed in [l|, in order to determine how it 
is affected by the new critical experimental information. 
Substantial changes with respect to the results presented 
in [l[ are expected. In fact, due to the strong anti- 
correlation existing among the two mixing angles (9 13 and 
014, the clear preference now emerged for a non-zero value 
of one of the two parameters (#13) should drastically re- 
duce the likelihood of the other one (#14) to be different 
from zero. Quantifying such a qualitative expectation 
appears particularly urgent in view of the numerous on- 
going projects of new experimental setups aimed at test- 
ing potential oscillations into sterile neutrinos (see for 
example [ij-|(|). 

The new landscape brings us to adopt a more conser- 
vative approach with respect to that espoused in [l| , as 
here our prime aim is to establish a robust estimate of #14 
independent of any assumption on the determinations of 
the reactor antineutrino fluxes. Indeed, their recent re- 
calculations [22[ , indicating an upward shift of about 
3% with respect to previous estimates, have not only re- 
fueled the interest around sterile neutrinos, but have also 
engendered an intense debate around possible systematic 
uncertainties, being common opinion that these may not 
be entirely under control. Having this issue in mind, we 
will treat the reactor data in a special way, minimizing 
the impact of the systematic uncertainties affecting the 
antineutrino fluxes. More specifically, in both the short- 
baseline reactor experiments (CHOOZ and D-CHOOZ) 
and the long-baseline ones (KamLAND), we will ignore 
the (flux dependent) total rate information, considering 
only that one provided by the energy spectral shape. 

This stratagem, although slightly limiting the con- 
straining power of the analysis, will render its results 
particularly robust. In fact, as discussed in [l[, the Kam- 
LAND analysis is quite sensitive to the reactor flux nor- 
malization. In particular, the indication in favor of non- 
zero 6*i3 (or 6*14) arising from the solar sector fluctuates 
between 1.3a and 1.8a, adopting, respectively, the old or 
the new (higher) fluxes Q. As a rule of thumb, we have 
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FIG. 1: Left panel: regions allowed after marginalization of 
the solar (Am 2 ot ,6i2) and atmospheric (Am 2 .^, #23) mass- 
mixing parameters by the solar sector data (diagonal bands) 
and LBL accelerator data (vertical bands). Right panel: re- 
gions allowed by their combination. The contours refer to 
Ax 2 = 1 (dotted line) and Ax 2 = 4 (solid line). 



verified that an upward (downward) 1% shift of the re- 
actor fluxes corresponds to a 0.1 5c increase (decrease) in 
the statistical significance of the preference for a non-zero 
electron neutrino mixing with ^3 (or 1/4). By removing 
the KamLAND total rate information from the analysis, 
we eliminate any dependency on the reactor flux nor- 
malization. In practice, with this procedure, the mixing 
angles 6*13 and #14 (and to a large extent also the "solar" 
mixing angle #12) are basically constrained by the solar 
data augmented 2 by the knowledge of the solar squared- 
mass difference Am 2 oi , whose high- precision determina- 
tion is preserved by retaining the KamLAND spectral 
shape information. 

Analogous considerations apply to the CHOOZ and 
D-CHOOZ experiments. Also in this case more (less) 
disappearance, and thus a preference of larger (smaller) 
values of #13 or 9n, is driven by higher (lower) reactor 
fluxes. Differently from KamLAND, however, the spec- 
tral information does not give any information on the 
relevant (atmospheric) mass-splitting AraJ tm , this being 
independently determined by the LBL — > disap- 
pearance searches performed at accelerators. It should be 
stressed that, in principle, the CHOOZ and D-CHOOZ 
spectral information could distinguish between the V3- 
driven (distorted) and 1/4-driven (undistorted 3 ) oscillated 



2 The solar data alone, without the "external" information on 
Am^ ( provided by the KamLAND spectral shape, would have a 
reduced sensitivity to all mixing angles. On the other hand, the 
KamLAND spectral shape provides little information on these 
last ones. 

3 It must be stressed that at the far detector (the only one cur- 
rently operational at the D-CHOOZ site) the oscillations driven 
by the new mass-mixing parameters (A m 2 , , 6 14) get com- 
pletely averag ed if Am 2 eu , > 0.1 eV 2 (see [23, HI). Therefore, 
in the region of the parameter space of current interest (con- 
fined to values of Am 2 eu) ~ 1 eV 2 ), we can safely assume that 



spectra, but its impact is negligible in practice since the 
expected distortions are very small (see the "Analysis C" 
in [HI). Indeed, the observation of such spectral distor- 
tions will be a challenge even for the next-generation of 
reactor experiments equipped with near detectors [25| . 
The achievement of this goal appears now even more im- 
portant in light of the opportunity of testing and distin- 
guishing standard and non-standard physics. 

Concerning the data sensitive to Am^, our analysis 
includes all the relevant solar and KamLAND data as de- 
scribed in detail in [l[ , but here the KamLAND absolute 
normalization is treated as a free parameter. As in [1] we 
have made the assumption that the additional mixing an- 
gles involving sterile neutrinos are null (#24 = #34 = 0). 4 

The regions allowed by the combined solar and Kam- 
LAND data represented by the diagonal bands in the 
left panel of Fig. 1 show no preference for non-zero mix- 
ing. This behavior, which is slightly different respect 
to that observed in [l[ (where we found a weak pref- 
erence for non-zero mixing), can be traced to the fol- 
lowing three factors: (I) The solar data taken alone give 
6»i3 = 6»i4 = as their best fit point; 5 (II) The KamLAND 
spectral shape taken alone does not show any preference 
for non-zero #13 or #14 ; 6 (III) The well-known interplay 
of KamLAND and solar data in pushing the #13 (0u) 
estimate upwards (see [H, @-H, HJ ) , so as to reduce the 
mismatch existing at #13 = 6*44 = among their (slightly 
different) determinations of the solar mixing angle 612, 
is now less effective since the KamLAND spectral shape 
has reduced sensitivity to this last parameter. 

Concerning the data sensitive to Am 2 tTW we incor- 
porate the LBL accelerator results as in account- 
ing for the —y disappearance searches performed 
at K2K [H and MINOS [33[, and the latest Va -> v e 
appearance results from MINOS Q and T2K gf. This 



the (Am 2 ew , 6i4)-induced oscillations are completely averaged 
with a consequent undistorted energy spectrum. Of course, the 
situation would be different at a detector located near to the 
reactor core (not operational at present), where non- negligible 
(Am 2 ew , 0i4)-induced spectral distortions are expected (see the 
discussion in [23j | ) . Finally, we remark that in the solar sector the 
new oscillations get averaged provided that Am 2 ellJ ^> Am 2 oi , 
as we have shown in 0. 

4 The assumption 24 = 6*34 = 0, implying in our parameterization 
(see [J) 11^4 = 0, is justified by the negative results of the short- 
distance disappearance searches performed in the Z/V — > chan- 
nel [26l . [27| , by the atmospheric data analyses [2g| | , and by the 
neutral current interaction searches performed by MINOS l29t . 
These last on es p rovide the stringent upper bound 024 < 7° at 
the 90% C.L. 129(1 . For such small values the 4i/-oscillation effects 
induced in LBL experiments, being (doubly) suppressed by the 
product I Ce4 1 1 UfjA I , would have a negligible impact in our analy- 
sis. In passing, we notice that it is for the same reason that the 
excess of the electron-like events observed in T2K and MINOS 
is not imputable to oscillations into sterile states. 

5 This feature has been also reported in other analyses @, @1 f° r 
what concern #13. 

6 Within a three-flavor framework the same behavior has been ob- 
served also in [3(ill , 
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dataset is insensitive to #14 and delimits the vertical band 
in the left panel of Fig. 1. To understand this point one 
should observe that z/4-driven —¥ v e appearance ef- 
fects are proportional to the mixing matrix element U^, 
which is set to zero in our analysis. 

The superposition (left panel of Fig. 1) of the two 
datasets sensitive respectively to Ato 2 o/ and Am 2 tm 
clearly evidences their complementarity in constraining 
the two mixing angles. Their synergy manifests quanti- 
tatively in their combination displayed in the right panel 
of Fig. 1. This provides the strong upper bound 



sin 2 6>i4 < 0.04 (90%C.L.) ; 



(1) 



which constitutes the main result of this report. For 
the sake of completeness, we mention that if we had in- 
cluded the total rate information from the reactor exper- 
iments we would have obtained a slightly weaker upper 
bound. For example, adopting the new (higher) fluxes' 
estimates the limit would become sin 2 6*14 < 0.05 at 90% 
C.L. In any case, the anti-correlation existing among the 
two mixing angles, characteristic of the solar sector, com- 
bined with the independent preference for non-zero #13, 
leads to a strong upper bound on #14, also destroying any 
weak preference for a non-zero value of this parameter. 
As an additional check of the robustness of the bound in 
Eq. (fT]) wc have verified that it is practically insensitive to 
the particular choice of the solar model used for the cal- 



culations. This is important in light of the yet unresolved 
"metallicity issue" and its connection with solar neutrino 
fluxes estimates (see [34J for an updated discussion of the 
topic) . 

We observe that the bound in Eq. ([I} is not incompat- 
ible with the estimates arising from the reactor (35| and 
gallium calibration [36H37I a nomalies. Rather, lying near 
their combined best fit [35|, it tends to select the lower 
part of the interval identified by such data. Probing such 
relatively low values of #14 with good precision should 
be the goal of any well-conceived experiment devoted to 
sterile oscillation searches. 

Finally, wc note that our limit is competitive with 
that one recently established in [38| using KARMEN and 
LSND i/ e -carbon cross sections, presenting the additional 
advantage of being independent of the new mass-squared 
splitting. This is a unique feature of the solar and reactor 
setups herein considered, where the new oscillations get 
completely averaged. 
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